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1.  Summary 


The  U.S.  Army  Research  Laboratory  (ARL),  Vehicle  Technology  Directorate  (VTD),  has 
recently  commissioned  a  state-of-the-art  Constant  Pressure  Flow  (CPF)  combustion  vessel 
facility  enabling  fundamental  spray  and  combustion  studies  of  critical  importance  to  the  mission. 
The  facility  operates  using  several  fuels  including  diesel,  JP-8,  Hydrotreated  Renewable  Jet 
(HRJ),  and  Iso-Paraffinic  Kerosene  (IPK)  and  can  carry  out  experiments  at  maximum 
temperature  and  pressure  conditions  of  1000  K  and  150  bar  while  maintaining  steady  conditions. 
Common  rail  and  Hydraulically-actuated  Electronically-controlled  Unit  Injector  (HEUI) 
injectors  provide  liquid  fuel  sprays  at  diesel  operating  conditions  typically  from  200  to  1800  bar. 
Spray  measurements  are  obtained  using  advanced  optical  diagnostics  such  as  Mie  scattering, 
shadowgraph,  and  schlieren  imaging.  Hence  the  facility  provides  a  critical  database  promoting 
technology  development  for  unmanned  aerial  vehicle  (UAV)  systems  (i.e..  Shadow,  Hunter,  and 
Gray  Eagle)  and  improvements  in  several  existing  power  focus  areas  (i.e.,  compact  components, 
increased  power  density  and  fuel  efficiency). 

In  response  to  this,  a  multiphase  flow  solver  CONVERGE,  developed  by  Convergent  Science, 
Inc.  (CSI),  has  been  adopted  to  enable  SD-Computational  Eluid  Dynamics  (CFD)  spray 
simulations  at  diesel  engine  operating  conditions.  The  solver  is  based  on  the  classical  Eularian- 
Eagrangian  formulation  and  blob-injection  wave  models  to  treat  the  spray  as  two  interacting 
fluids.  Over  the  years,  conducting  spray  simulations  have  demonstrated  to  be  of  paramount 
importance  in  the  research  process.  Some  important  features  include  its  suitability  in  conducting 
extensive  cost-effective  parametric  and  optimization  studies,  its  ability  to  provide  a  complete 
history  of  multidimensional  resolved  information  for  every  variable,  and  its  ability  to  artificially 
separate  specific  subprocesses  that  would  inherently  interact  in  an  experiment.  Although  CED- 
based  tools  have  steadily  improved  and  continue  to  provide  higher  quality  models,  a  vast  number 
of  experiments  will  remain  necessary  to  guide  and  control  simulation  results,  to  validate  and 
calibrate  existing  models,  and  in  driving  the  development  of  next-generation  models. 
Consequently,  a  much  better  analysis  for  the  interpretation  of  complex  results  will  be  available, 
if  both  numerical  and  experimental  studies  are  conducted  in  parallel. 

Historically,  the  CED  community  has  benefited  from  spray  measurements  that  have  driven  the 
model  development  process.  Validation  studies  will  reveal  the  suitability  of  the  modeling 
assumptions  (physical  models),  stability  of  the  spatio-temporal  numerical  technique  (numerical 
methods),  and  calibration  of  model  parameters  (turbulence,  breakup,  combustion  constants)  that 
is  required  to  carefully  benchmark  a  simulation.  In  this  report,  several  studies  are  conducted  at 
diesel  engine  operating  conditions  to  demonstrate  the  modeling  capability  of  CONVERGE.  The 
databases  used  are  part  of  the  Engine  Combustion  Network  (ECN)  library.  The  results 
demonstrate  the  suitability  of  the  models  used,  physical  assumptions,  and  boundary  conditions. 
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Spray  simulations  are  in  good  agreements  with  experiments,  providing  eonfidenee  in  modeling 
proeedure,  and  enabling  further  studies. 


2.  Introduction 


In  many  combustion  systems,  fuel  atomization  and  the  spray  breakup  process  play  an  important 
role  in  determining  combustion  characteristics  and  emissions  formation.  Due  to  the  ever  rising 
need  for  better  fuel  efficiency  and  lower  emissions,  the  development  of  a  fundamental 
understanding  of  its  process  is  essential  and  remains  a  challenging  task.  In  the  U.S.  military  air 
and  ground  vehicle  applications,  the  characterization  of  engines  operating  with  JP-8  fuels  is 
critical  in  developing  new  vehicle  propulsion  technologies  and  enhancing  existing  compact  fuel- 
to-electric  power  systems.  Hence,  the  control  of  spray  characteristics  will  remain  of  critical 
priority  in  enabling  the  Army  mission  of  driving  new  technologies,  optimizing  engine 
performance,  while  facilitating  fuel  logistics  operations. 

In  modem  engines,  the  spray  process  begins  with  a  high-pressure  liquid  fuel  injected  into  a 
vessel  at  high  temperature  and  pressure  conditions.  A  liquid  core-jet  is  formed;  its  penetration 
length  is  influenced  by  several  factors  including  nozzle  orifice,  internal  nozzle  conditions,  fuel- 
air  turbulent  mixing,  external  ambient  conditions,  and  physical  properties.  The  interaction  of  the 
fuel  jet  with  its  surroundings  promotes  liquid  atomization  into  fine  droplets  as  a  result  of 
aerodynamic  instabilities  and  surface  tension  forces.  The  continuation  of  the  process  leads  to 
further  breakup,  droplet-droplet  interactions,  droplet-wall  interactions,  and  evaporation.  If 
critical  conditions  are  present,  the  fuel  may  self-ignite  leading  to  a  combustion  event  and 
generating  a  specific  amount  of  power  output.  Research  studies  providing  access  to  spray 
measurements  will  remain  critical  in  enabling  new  understanding  of  the  spray  process. 

Remarkable  progress  has  been  made  characterizing  the  spray  behavior  in  both  experimental  and 
numerical  studies.  Experimentally,  Siebers  (1)  has  studied  the  behavior  of  evaporating  diesel 
sprays  and  has  reported  the  effects  of  several  operating  condition  on  the  liquid  and  vapor  spray 
penetrations,  kinematics,  and  mixture  fraction.  It  was  reported  that  liquid  length  decreases 
linearly  with  injector  orifice  diameter,  temperature,  and  increases  with  fuel  volatility  or 
temperature,  and  that  it  is  weakly  dependent  on  injector  pressure.  In  a  cross-institutional  study. 
Picket  et  al.  (2)  has  presented  the  operational  details  of  lab  instmmentation,  and  experimental 
conditions  across  several  platforms  in  an  effort  to  enhance  the  fidelity  of  spray  diagnostics. 

Payri  et  al.  (5)  has  studied  the  fuel-temperature  effect  in  nonreacting  and  reacting  diesel  sprays 
using  a  novel  injector,  and  imaging  diagnostics  for  liquid  phase  penetration,  light-off  length,  and 
ignition  delay  measurements.  Payri  et  al.  (5)  reported  a  clear  influence  of  the  injector  body 
temperature,  real  fuel  temperature,  on  nonreactive  and  in  a  lesser  degree  to  reacting  sprays.  More 
recent  studies  have  focused  on  resolving  the  structure  of  the  spray’s  liquid  core  by  use  of  x-ray 
techniques  and  medical  imaging  {4-6).  Wang  et  al.  {4)  used  x-ray  phase  contrast  imaging  to 
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study  the  near  nozzle  atomization  proeess  of  air  assisted  water  sprays.  Key  observations  were 
made  about  the  atomization  proeess  at  high-We  numbers,  sueh  as  jet  narrowing,  spray  breakup, 
and  the  tracking  of  the  mass  volume  fraction.  Coletti  et  al.  (6)  has  introduced  x-ray  computed 
tomography  (CT)  technique  to  provide  detailed  information  of  the  spray  dense  region.  Medical 
software  is  used  for  imaging  to  reconstruct  a  three-dimensional  time-averaged  distribution  of 
liquid  mass  fraction  from  the  intact  core  to  the  dilute  spray  region.  This  progress  has  provided 
new  insights  into  the  spray  including  the  near  nozzle  region  improving  our  understanding  and 
driving  the  generation  of  enhanced  spray  models. 

Over  the  years,  3D-CFD  multiphase  solvers  have  been  developed  largely  based  on  Lagrangian 
particle  tracking  methods  coupled  with  wave  type  spray  models  and  appropriate  source  terms.  In 
a  study  of  spray  induced  turbulence,  Banerjee  and  Rutland  (7)  has  used  a  traditional  code 
(KIVA)  with  a  selection  of  modem  turbulence  models  on  various  grids.  The  study  reported  on 
the  resolution  study  conducted  with  large-eddy  simulation  (LES)  methods  leading  to  CFD  mesh 
criteria  based  on  the  breakdown  of  large  flow  coherent  stmctures.  Also  using  the  KIVA  solver, 
Bharadwaj  et  al.  (5)  presented  LES  results  for  a  particle  laden  gas  jet  and  nonevaporative  diesel 
sprays  obtaining  good  agreements  with  Argonne  National  Laboratory  (ANL)  x-ray  measured 
data.  The  objective  was  to  demonstrate  the  capability  of  LES  in  predicting  two-phase  interaction 
while  modeling  of  the  subgrid  scale  droplet  effects  on  turbulent  kinetic  energy.  More  recently, 
Senecal  et  al.  (P)  has  demonstrated  grid-convergent  properties  using  the  state-of-the-art 
CONVERGE  solver  and  provided  guidelines  to  benchmark  against  ECN  spray  measurements.  In 
an  effort  to  capture  injector  cavitation  and  turbulence  effects  on  spray  atomization  models,  Som 
et  al.  {10)  introduced  the  Kelvin-Helmholtz -Aerodynamics  Cavitation  Turbulence  (KH-ACT) 
model  and  demonstrated  its  fidelity  through  validation  across  several  spray  configurations.  It  was 
observed  that  the  inclusion  of  cavitation  and  turbulence  enhances  primary  atomization 
decreasing  spray  droplet  sizes,  reducing  penetration  lengths,  and  increasing  radial  dispersion.  As 
a  consequence,  the  model  provides  better  agreements  with  measurements  of  non-evaporating  and 
evaporating  sprays.  In  2013,  Senecal  et  al.  {11)  extended  his  previous  spray  study  to  transient 
LES  methods,  recalibration  of  spray  parameters,  and  the  effect  of  cycle  to  cycle  averaging  on 
spray  diagnostics.  A  similar  study  was  conducted  by  Xue  et  al.  {12)  focusing  on  grid- 
convergence  studies  with  comparison  to  ECN  data  denoted  as  Spray-A  and  Spray-H  and  hence 
developing  grid  criteria.  These  studies  have  enabled  practical  engineering  simulations  leading  to 
enhanced  understanding  and  advancements  in  technology  developments  and  research  process. 

In  this  report,  several  3D-CFD  validation  studies  are  conducted  to  demonstrate  the  modeling 
capability  of  the  CONVERGE  solver.  Two  databases  were  selected  providing  a  wide  range  of 
operating  condition  relevant  to  diesel  engine  configurations.  The  databases  used  are  from  the 
ECN  and  are  representative  of  present  research  studies  providing  a  synergistic  approach  to 
research.  The  report  is  organized  in  the  following  manner.  A  brief  literature  review  on 
experimental  and  numerical  spray  studies  has  been  presented  with  a  focus  on  recent  trends.  A 
description  of  the  governing  equations,  time  advancement  method,  turbulence,  and  spray  model 
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is  presented  to  establish  the  aecuraey  of  the  solver  and  its  numerical  methodology.  The  results 
are  then  discussed  in  two  sections  where  comparisons  are  made  to  the  ECN  measurements  at 
various  conditions.  The  ECN  databases  used  correspond  to  the  classical  Spray  A  conditions  from 
Sandia  National  Eaboratories,  and  CMT-Motores  study  of  fuel-temperature  effects  on  spray 
penetration  behavior.  Conclusions  and  future  efforts  are  discussed. 


3.  Numerical  Method 


The  CONVERGE  3D-CED  solver,  developed  by  CSI,  has  been  adopted  in  this  study  to  perform 
detailed  spray  simulations  at  realistic  engine  operating  conditions.  CONVERGE  is  a 
compressible  Navier  Stokes  solver,  it  is  based  on  a  first  order  predictor-corrector  (Pressure 
Implicit  with  Splitting  of  Operator  [PISO])  time  integration  scheme,  and  a  choice  of  second  or 
higher  order  finite  volume  schemes  for  spatial  discretization.  It  features  a  nonstaggered, 
collocated,  computation  grid  framework  using  a  Rhie-Chow  interpolation  technique  to  avoid 
spurious  oscillations.  An  efficient  geometric  multigrid  treatment  is  used  to  solve  the  pressure 
equation,  and  a  parallel  computing  implementation  is  based  on  implementations  of  either 
OpenMP  or  Message  Passing  Interface  (MPI)  protocols.  It  provides  the  option  of  increasing 
resolution  locally  through  static  fixed-grid  embedding  and  dynamically  through  Adaptive  Mesh 
Refinement  activated  through  user  specified  criteria.  Additionally,  it  uses  state-of-the-art 
Eularian-Eagrangian  spray  models,  and  a  parallel  detailed  chemistry  solver  for  combustion  that 
can  be  fully  coupled  with  CHEMKIN  databases.  CONVERGE  also  provides  a  choice  between  a 
number  of  modeling  options  for  the  treatment  of  turbulence,  including  direct  numerical 
simulation  (DNS),  EES,  and  Reynolds-averaged  Navier-Stokes  (RANS). 

3,1  Governing  Equations 

In  this  study,  the  gas  phase  is  described  using  the  Eavre-averaged  Navier-Stokes  Equations,  and 
the  RNG  k  —  6  model  for  low-resolution  turbulence  studies  has  been  adopted.  The  compressible 
system  of  transport  equations  for  mass,  momentum,  energy,  and  species  transport  are  presented 
here  in  a  RANS  framework. 


Conservation  of  Mass: 


^  _l_  dpuj  _  Q 
dt  dxj 


Conservation  of  Momentum: 


dpui  dpUiUj 
dt  dXj 


dP  9  r  PdUj  dUj\ 
dxi  dxi  ^  \9xj  dxi) 


(1) 

(2) 
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The  last  term  is  denoted  as  the  Reynolds  stresses  of  the  system,  Tjy  =  —pu[Uj.  It  needs  to  be 
modeled  to  provide  mathematieal  elosure  and  to  aeeount  for  turbulenee  effeets.  Turbulent 
viscosity  formulation  in  RANS  typically  require  the  modeled  Reynolds  stress  for  the  Standard 
k  —  6  and  RNG  to  be  represented  as, 


'i] 


=  pu[u'j  =  ZptSij  -^Sij  (p/c  + 


(3) 


Hence,  turbulent  viscosity  is  defined  as  pt  —  C^p^k^je),  where  k  and  e  are  calculated  classically 
by  transporting  two  representative  equations.  Additionally,  there  is  a  model  constant  that 
typically  carries  the  value  of  =  0.09,  and  Sij  represents  the  mean  strain  rate  tensor.  The 
transport  equations  for  the  k  and  e  RANS  model  are  given  in  order  as. 


dpk 

dt 


dpUik 

dxi 


dui 


(4) 


and. 


^  +  £(0£)  =  A ('jLil')  _  c,,pe ^  +  f C,3 

dt  dXi  dXj\PrgdXjJ  dxi  \  dxj 


(5) 


where  the  definition  of  R  will  vary  depending  on  the  choice  of  RANS  turbulence  model.  It  is 
defined  as  R  =  0  for  the  Standard  k  —  e  model  and  R  =  —  r\lr\p^e^  jiX.  —  pri^)k  for 

RNG  k  —  6  model. 

Note  that  the  source  terms,  Sg  ,  are  included  due  to  the  Lagrangian-Eularian  formulation  and  are 
important  in  describing  the  particle  interaction  with  its  surrounding.  It  is  defined  as. 


where  Np  is  the  number  of  drops  in  a  parcel,  V  is  the  cell  volume,  ii[  is  the  gas  phase  turbulence 
fluctuation,  and  the  summation  is  carried  over  all  the  parcels  in  each  cell.  Note  that  the  drag 
force  on  each  drop  is  defined  as,  =  Fdrag.i'^il'^i  + 

Heat  and  species  mass  transfer  is  treated  using  the  compressible  form  of  the  energy  equation  and 
the  species  mass  transfer  equations.  Note  that  the  species  transport  solves  for  the  mass  fraction, 
Tfc,  of  all  species  in  the  specified  system. 


Conservation  of  Energy: 

dpe  ^  diTjpe 
dt  dxj 


PF  Yik 


dxj) 


+  s 


(V) 


5 


Conservation  of  Species: 


dpk  _|_  ^Pk'^j 
dt  dxj 


(8) 


where  e  is  the  specific  internal  energy,  K  is  the  conductivity,  is  the  species  enthalpy,  D  is  the 
mass  diffusion  coefficient,  and  the  source  term.  For  turbulent  cases,  the  mass  diffusion 
coefficient  is  calculated  by  Df-  —  Vf-/SCf-  in  terms  of  turbulent  viscosity  and  Schmidt  numbers. 


3.2  Time  Advancement 


Equations  1  and  2  are  solved  numerically  using  a  classical  predictor-corrector  scheme  in  which 
the  velocity  field  is  first  integrated  using  the  Navier  Stokes  equation  2  and  then  corrected  to 
enforce  mass  conservation  (equation  1)  using  a  modified  pressure.  This  is  achieved  using  the 
modified  PISO  algorithm  first  introduced  by  Issa  et  al.  (13). 

An  illustration  of  the  PISO  predictor  corrector  sequence  is  presented  below  (note  that 
H*  represents  convection,  diffusion,  S  accounts  for  implicit/explicit  sources,  and  *  denotes  the 
intermediate  steps)  for  a  one-time  correction  procedure. 


Predictor  step: 


Corrector  steps: 


Update  step: 


dt 


_ 

dt 


♦  ** 

P  Uj 

dt 


dt  dxi  ^ 

(p*_pn)0 


dxidxi 


(p*  _  pn)  _  ^  (dp^u*  _\J_ 

dt^  V  dXi  )  dt 


*  ** 

P  Uj 

dt 


P'^u*  _ 
dt 


dXi  ^ 


(9) 


(10) 

(11) 


(12) 


where  dt  is  the  numerical  integration  time  step.  Once  the  pressure  is  solved,  the  velocity  is 
recalculated  according  to  the  updated  momentum  equation. 
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3,3  Spray  Modeling 

Spray  modeling  is  treated  using  the  “blob”  injection  method  of  Reitz  and  Diwakar  {14).  Blobs  of 
a  characteristic  size  are  injected  following  a  statistical  distribution  into  the  computational 
domain.  Primary  and  secondary  breakups  are  subsequently  simulated  based  on  the  Kelvin- 
Helmholts  (KH)  and  Rayleigh-Taylor  (RT)  instability  methods.  Note  that  the  breakup  length  is 
not  determined  a  priori  (breakup  length  concept)  and  is  calculated  as  a  part  of  the  solution. 

In  the  KH  wave  model,  atomization  is  treated  using  stability  analysis  for  liquid  jets.  The  breakup 
of  the  injected  blobs  and  resulting  drops  of  radius  is  calculated  by  assuming  that  the  drop 
radius  is  proportional  to  the  wavelength  of  the  fastest  growing  unstable  surface  wave  . 

It  is  written  as, 


r  —  BoAi^fq 


(13) 


where  Bq  is  a  model  constant.  The  droplet  size,  and  its  change  of  radius  is  in  the  following  way. 


^  _  {rp-r) 

dt  TxH 


(14) 


where  the  breakup  time  constant,  ,  is  calculated  as, 

3. 7265^7^0 


(15) 


and  the  maximum  growth  rates  Hkh  and  corresponding  wavelengths  have  been  simplified 
and  defined  as  follows. 


fpia^\  _  0.34+0.38Vl/e^-^ 
\  a  J  ~  (1+Z)(1+1.4T0'6) 


(16) 


and, 

^  ^  g  Q2  (l+0-45Z°-^)(l+0-4r°  ") 

a  '  (l+0.87Wej®^)°'® 


(17) 


where, 

Z  =  We^-^/Rei ,  T  =  ^  Wei  =  PiU^ala  ,  WCg  =  PgU^aja,  and  Rei  =  Ua/vi 


The  present  RT  mechanism  formulation  includes  viscosity  variations  in  the  growth  rate  equation, 

(18) 


0)^7 


\Pl+Pg) 


\Pl+Pg)  Pl+Pg  \Pl+PgJ 


where  is  the  wavenumber,  pi  is  the  liquid  viscosity,  Pg  is  the  gas  viscosity,  pi  is  the  liquid 
density,  Pg  the  gas  density,  a  is  the  deceleration  of  the  drop,  and  cr  is  the  liquid  surface  tension. 
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The  wave  number  corresponding  to  the  maximum  growth  rate  Knj  =  2n/Aiij  is  calculated 
through  a  bisection  method  with  equation  18. 

The  value  is  updated  to  calculate  the  maximum  growth  rate  The  predicted  RT  model  drop 
size  is  then  expressed  as, 

Vrt  =  CrjArr  (19) 

where  Crj  is  the  model  constant,  and  Arr  is  the  predicted  RT  wavelength.  In  summary  the  spray 
model  constants  used  in  this  study  are  defined  as:  Sq  —  0-6 , 5i  =  12  ,  Crj  —  0.1 ,  and  —  1. 

Lastly,  other  spray  processes  that  were  modeled  include  droplet  distortion  and  drag,  droplet 
interactions  in  terms  of  collision  and  coalescence,  turbulent  dispersion,  and  evaporation.  The  two 
phases  are  coupled  through  the  exchange  of  mass,  momentum,  and  energy,  represented  in  the 
appropriate  source  terms  in  the  gas-phase  equations.  These  models  are  described  in  detail  in  a 
recent  review  by  Bravo  and  Kweon  (15)  and  the  reader  is  referred  to  that  article. 


4.  Results  and  Discussions 


Measurements  selected  from  two  databases  were  utilized  to  benchmark  the  CFD  solution  at 
various  conditions.  The  databases  can  be  readily  extracted  from  the  ECN  Internet  library 
operated  by  Sandia  National  Laboratories,  Combustion  Research  Facility  and  available  at 
http://www.sandia.gov/ecn/.  The  network  provides  access  to  detailed  spray  measurements 
conducted  through  the  use  of  multi-institutional  combustion  vessels  (CPF  or  Constant  Volume 
Prebum  [CVP]  chambers)  at  various  ambient  gas  compositions.  Some  of  the  ECN  affiliates 
across  several  institutions  that  have  conducted  spray  studies  include:  Sandia  National 
Laboratories  (U.S.),  lEP  (Erance),  Universidad  Politecnica  de  Valencia  (Spain),  Caterpillar 
(U.S.),  and  Eidenhoven  University  of  Technology  (Netherlands). Their  efforts  have  contributed  to 
the  breadth  of  the  library  by  conducting  reacting/non-reacting  spray  studies  under  several 
operating  conditions,  and  fuels  (i.e.,  gasoline,  diesel) 

Recently,  several  efforts  have  been  reported  attempting  to  standardize  experimental 
configurations  and  procedures.  An  ideal  configuration  was  created,  denoted  Spray  A,  where  the 
influence  of  cross  institution  facilities  and  conditions  was  reported.  This  includes  the  influence  of 
injectors,  combustion  vessels,  and  optical  diagnostic  tools.  Hence,  Spray  A  is  a  well-documented 
configuration  and  has  been  selected  as  the  baseline  case  study. 
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4,1  ECN,  Sandia  Spray  A  Conditions 

In  this  configuration,  the  experiment  is  carried  out  in  a  CVP  vessel  chamber  faeility  at  Sandia 
National  Laboratories  as  shown  in  figure  1 .  The  desired  pressure  and  temperature  conditions  are 
aehieved  by  spark  igniting  a  premixed  gas  mixture  to  eompletion.  The  eombustion  produets  eool 
over  a  period  of  700-2500  ms  due  to  interaetion  with  the  vessel  walls,  with  the  ehamber  pressure 
deereasing  slowly.  Once  the  target  eonditions  are  deteeted,  the  fuel  injector  is  triggered  and  the 
spray  process  begins. 


Figure  1.  Sandia  National  Laboratories  (a)  CVP  facility,  (b)  schematic  of  CVP  setup. 


The  experimental  faeility  also  provides  full  optical  access  for  line  of  sight  or  orthogonal  optieal 
diagnosties  as  seen  in  figure  la.  Figure  lb  shows  the  setup  with  the  positioning  of  the  high- 
pressure  eommon  rail  fuel  injector,  number  of  spark  plugs,  and  fan  location.  To  prevent  wall 
impingement  effeets  interfering  during  spray  diagnostic  process,  the  characteristie  length  of  the 
vessel  is  typically  designed  to  be  100  mm  on  eaeh  side.  In  order  to  prevent  condensation  of 
combustion  products  on  the  windows  the  vessel  walls  are  heated  eleetrieally  to  engine  surface 
temperatures.  Multiple  spark  plugs  are  used  to  provide  eonsistent  ignition  of  a  prebum  lean 
mixture.  Further  details  of  the  vessel  geometry  and  diagnosties  are  presented  in  the  review  by 
Meijer  et  al.  {16). 

Spray  A  conditions  are  provided  in  table  1 .  The  physical  description  corresponds  to  an 
evaporating  fuel  spray  with  0%  oxygen  content  (nonreacting),  developing  at  diesel  engine 
ambient  eonditions.  A  single  hole,  modern  eommon  rail  injeetor  with  an  injector  diameter  of 
90  p  (Boseh  CRIN  2.4)  is  used  at  typical  diesel  injection  pressures.  A  single  eomponent  diesel 
surrogate  fuel  (i.e.,  n-dodeeane)  is  used  due  to  its  extensively  well-characterized  ehemieal  and 
physical  properties.  Detailed  and  redueed  mechanisms  for  n-dodecane  are  readily  available  from 
the  literature  (77). 
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Table  1.  Injector  and  spray  conditions  (Spray  A). 


Sandia  Conditions 

Valne 

Sandia  Conditions 

Valne 

Fuel 

n-dodecane 

Injection  Pressure  (Mpa) 

150 

Ambient  composition 

0%  Oxygen 

Fuel  Temperature  (K) 

363 

Ambient  temperature  (K) 

900 

Nozzle  Diameter  (mm) 

0.09 

Ambient  density  (kg/m^) 

22.8 

Injection  Duration  (ms) 

1.5 

Number  of  injector  holes 

1  (axial) 

Injection  mass  (mg) 

3.5 

Numerical  modeling  with  CONVERGE  was  conducted  by  using  the  multiphase  Lagrangian 
particle  tracking  methods  with  wave-based  spray  submodels  and  the  RANS  methodology  as 
described  in  section  3.  There  are  several  important  procedures  to  be  mindful  of  when  running 
spray  simulations.  An  Adaptive  Mesh  Refinement  (AMR)  procedure  must  be  adopted  to  control 
the  refinement  levels  in  the  vicinity  of  the  spray.  In  this  study,  AMR  is  activated  for  the  velocity 
field  and  its  use  is  critical  in  keeping  the  cell  counts  at  realistic  operating  levels. 

Note  also  that  fixed-grid  embedding  (in  spray  liquid  core  region)  is  used  to  define  the  reference 
minimum  cell  sizes,  dx,  and  the  embedding  scale.  The  parameters  are  related  through  the  base 
grid  cell  size,  dxij^se  ?  follows:  dx  =  dxij^se  ^  2 /{embed  scale^).  The  coarsest  base  grid 
cell  size  of  2  mm  and  an  embed  scale  of  3  corresponds  to  dx  —  0.25  mm. 

When  conducting  grid  refinements,  a  fixed  number  of  injected  particles  will  lead  to  an  artificial 
increase  of  the  liquid  penetration  length  (P).  This  is  due  to  the  decrease  of  relative  gas  to  liquid 
mass  that  occurs  in  each  cell  as  the  mesh  is  refined.  Overlooking  this  can  artificially  cause  the 
gas  velocity  to  approach  the  particle  velocity  at  every  cell,  decreasing  particle  drag.  Table  2 
documents  the  number  of  resolution  levels,  denoted  as  embedding  scale,  and  the  critical  number 
of  injected  particles  to  achieve  converged  results.  Note  that  to  avoid  going  below  a  critical 
injected  parcel  number  of  3000,  the  1  mm  and  2  injected  parcel  cases  are  identical. 


Table  2.  CFD  embed  scale  and  number  of  injected  particles. 


Embed  Scale 

dx 

(mm) 

Injected  Parcels 

0 

2 

3125 

1 

1 

3125 

2 

0.5 

12,500 

3 

0.25 

50,000 

4 

0.125 

200,000 

5 

0.0625 

800,000 

Eigure  2  shows  the  model  convergent  properties  by  inspection  of  the  multiphase  penetration 
profiles.  Spatial  grid-refinement  studies  were  conducted  using  the  grid  and  Eagrangian  parcel 
parameters  listed  in  table  2.  In  figure  2a,  the  liquid  penetration  diagnostic  shows  good  agreement 
with  the  experimental  data  and  convergent  behavior.  Eor  the  coarse  case  of  dx  =  2  mm  and 
embed  scale  of  zero,  liquid  penetration  is  over-predicted  at  peak  by  a  factor  of  two.  No 


10 


noticeable  improvement  is  obtained  with  one  additional  level  of  grid  refinement.  However,  the 
behavior  drastically  improves  as  the  embedding  seale  elevates  from  a  level  of  three  to  five,  and 
finally  converges  to  the  experimental  data.  Figure  2b  shows  a  similar  trend  with  the  vapor 
penetration  length  with  the  solution  eonverging  as  the  grid  is  refined.  In  this  ease,  however,  only 
eonvergenee  is  obtained  while  the  model  solution  at  that  refinement  level  shows  discrepancies 
with  the  measured  data.  Figure  2c  shows  the  convergence  behavior  of  the  peak  liquid  penetration 
with  resolution  level  providing  a  guideline  for  grid  criteria.  It  is  important  to  note  that  due  to  the 
steady  nature  of  the  RANS  turbulenee  model  no  cyele-to-eycle  averaging  was  conducted  and  a 
single  shot  analysis  is  sufficient. 


Figure  2.  Spray  A  convergence  properties  at  various  embed  scales  corresponding  to  dx  =  2,  1,  0.5,  0.25,  0.125, 
0.0625  mm:  (a)  liquid  penetration  vs.  time,  (b)  vapor  penetration  vs.  time,  (c)  peak  liquid  penetration 
vs.  embed  scale. 
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In  figure  3,  an  embedding  seale  of  three  was  seleeted  beeause  it  provided  a  good  eompromise 
between  aeeuracy  and  eomputational  eost.  Also  shown  in  figure  3  are  multiphase  penetration 
eomparisons  with  Sandia,  CMT-Motores  (Spain),  and  IFF  (Franee)  measurements  for  additional 
referenee  (2).  The  liquid  penetration  profile  in  figure  3a  shows  two  elearly  marked  regions.  The 
regions  are  denoted  by  a  transient  region  oeeurring  up  to  a  time  of  0.25  ms  and  after  this 
threshold  a  steady  region.  Note  a  penetration  overshoot  within  the  transient  region  of  figure  3a. 
Although  refining  the  grid  does  attenuate  this  over-shoot  effeet  (see  figure  2a),  it  does  not 
eliminate  it.  This  effeet  is  thought  to  oeeur  due  to  the  present  RANS  turbulenee  treatment  and  its 
behavior  with  high-fidelity  LES  or  DNS  models  and  will  be  addressed  in  future  works. 


Figure  3.  Multi-institutional  experimental  comparison  of  Spray  A  penetration  profiles:  (a)  liquid  penetration, 
(b)  vapor  penetration.  Embed  scale  of  three,  dx  =  0.25  mm. 

Figure  3b  shows  the  vapor  penetration  profile  for  Spray  A  at  an  embedding  seale  of  three.  It  is 
noted  that  converged  profile  deviates  from  the  experimental  data,  as  was  also  the  case  in  figure 
2b.  This  effect  has  been  observed  previously  by  Senecal  et  al.  (P)  where  he  attributed  this 
discrepancy  to  current  uncertainties  in  measurement  techniques.  However  the  penetration 
velocity  and  rate  of  vapor  penetration  (slope  of  the  profiles)  is  in  close  agreement. 

Figure  4  shows  comparison  with  the  experimental  images  reported  by  Pickett  et  al.  (2)  for  the 
same  conditions.  The  measurements  were  obtained  using  Mie  scatter  optical  diagnostics  and  the 
images  are  time  averaged.  Note  that  the  reporting  discrepancies  with  the  IFP  data  (peak  in 
penetration  profile)  has  been  reported  and  thought  to  be  due  to  differences  in  laser  optical 
diagnostic  methodologies  (orientation  of  laser  light). 
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Figure  4.  Mie  scattering  imaging  sequence  at  IFF  and  Sandia  and  comparison  to  ARL’s  CFD  model  for 
Spray  A  conditions. 
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Figure  5a  and  5b  show  the  Eularian  gas  phase  velocity  and  fuel-mass  fraction  (n-dodecane) 
contour  variations  surrounded  by  the  vapor  mass  fraction  line,  Yf  —  0.01.  Also  shown  are  the 

Lagrangian  liquid  phase  particles  convecting  downstream  subject  to  a  uniform  heating  load  and 
shear  forces.  Note  that  peak  values  will  originate  from  the  nozzle  orifice  rapidly  decreasing  in 
magnitude  due  to  mixing  process  with  the  stagnant  environment.  Time  intervals  corresponding  to 
0.5,  1.0,  and  1.5  ms  are  selected  to  show  the  spray  process. 


0 


velocity  (m/s) 

50  100150  2^^0300^0 


yNC12H26 

0.00  0.06  0.12  0.18  0.24  0.30 


Figure  5.  Spray  A  visualization  (a)  contour  of  velocity  and  (b)  contours  of  fuel-mass  fraction.  Note  that 
spray  is  surrounded  by  the  vapor  line  at  =  0.01. 
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4,2  ECN,  CMT  Spray  A  conditions  (Fuel-Temperature  Effects) 

In  this  configuration,  the  experiment  is  earried  out  in  a  CPF  test  ehamber  that  has  the  capability 
to  reach  peak  conditions  of  1000  K  and  15  MPa,  figure  6a.  It  does  not  require  a  pre-bum  phase  to 
obtain  this  eondition  and  rather  it  is  eomprised  of  four  subsystems  (1)  gas  eompressor,  (2)  gas 
heater,  (3)  test  vessel,  and  (4)  control  system  to  achieve  nearly  quiescent  and  a  steady 
thermodynamic  conditions.  A  gas  initially  stored  in  a  high-pressure  reservoir,  flows  continuously 
through  the  vessel  via  volumetrie  eompressors.  Eleetrieal  heaters  are  used  to  inerease  the 
temperature  of  the  flowing  gas  up  to  the  seleeted  temperature.  A  control  system  (closed  loop 
proportional-integral-derivative  [PID]  controller)  adjusts  for  both  the  chamber  pressure  and  the 
power  of  the  heater  to  obtain  uniform  conditions  for  the  test.  To  minimize  temperature 
inhomogeneity  in  the  testing  region  the  thinner  inner  walls  are  covered  with  an  insulating  layer 
and  have  the  function  of  reducing  the  heating  of  the  outer  shell.  Note  the  initial  conditions  differ 
from  Sandia’s  combustion  vessel,  since  there  are  no  eombustion  products  present. 


Figure  6b  shows  the  optieal  setup  used  for  Mie  seattering  technique  with  frontal  illumination. 
The  spray  region  is  illuminated  with  a  eontinuous  Xenon-arc  lamp  where  the  seattered  light  is 
recorded  by  way  of  a  high-speed  CMOS  camera  (Photoron  Fastcam). 


(b) 


Figure  6.  CMT-Motores,  Universidad  Politecnica  de  Valencia:  (a)  CPF  facility,  (b)  schematic  of  CPF  setup. 

Baseline  Spray  A  conditions,  as  provided  in  table  1,  are  also  invoked  in  this  study.  It  is  important 
to  note  the  difference  in  setup  between  the  two  studies.  In  the  present  study,  there  is  a  strong 
emphasis  on  characterizing  fuel-temperature  effects,  and  as  such,  fuel-coolant  temperatures  are 
varied  accordingly.  In  addition,  the  vessel  gas  temperature  is  allowed  to  vary  from  820  to  900  K, 
as  well  as  the  injection  pressure  from  50,  100,  to  150  MPa.  Vessel  conditions  are  given  in  table 
3. 
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Table  3.  CMT  vessel  conditions. 


CMT  Conditions 

Value 

CMT  Conditions 

Value 

Fuel 

n-dodecane 

Injection  pressure  (Mpa) 

50,  100,  150 

Ambient  composition 

100%  Ny, 

Fuel  temperature  (K) 

363 

Ambient  temperature  (K) 

820,  900 

Nozzle  diameter  (mm) 

0.09 

Ambient  density  (kg/m^) 

22.8 

Injection  duration  (ms) 

4.5 

Number  of  injector  holes 

1  (axial) 

Injection  mass  (mg) 

3.5 

Note  that  there  is  a  strong  distinetion  from  the  fuel-eoolant  temperature  and  the  real  fuel 
temperature.  The  relationship  between  fuel  eoolant  and  real  temperature  has  been  studied 
through  the  use  of  a  dummy  injeetor  test  rig.  Beeause  there  is  no  fluid  flow  through  the  dummy 
injeetor  the  outer  nozzle  hole  has  been  shut,  while  the  injeetor  needle  eontains  an  axially  drilled 
hole  with  the  upstream  part  to  insert  a  thin  thermoeouple.  Temperature  readings  are  made  axially 
up  to  3 3 -mm  upstream  from  the  injeetor  sae  to  provide  temperature  readings  in  this  region  and 
study  the  relationship  between  both  lines.  For  further  referenee  on  the  dummy  injeetor  set  up  and 
instrumentation,  see  Payri  et  al.  (5). 

Table  4  presents  the  measured  fuel  temperature  at  the  injeetor  sae  volume  (Tf)  for  two  vessel 
ambient  eonditions,  T^mb  —  ^20  and  900  K.  Three  eoolant  eonditions  (T^)  were  tested  for  eaeh 
vessel  eondition  and  the  sae  fuel  temperature  was  determined  through  a  sliding  thermoeouple 
instrumentation.  In  this  study,  the  sae  fuel  temperature,  Tf  ,  has  been  used  as  the  best 
approximation  to  the  CFD  model  fuel-temperature  boundary  eondition. 


Table  4.  Coolant  and  injector  sac  fuel  temperature. 


Tamb  (K) 

Tamb  (K) 

Tc(K) 

Tamb  (K) 

Tamb  (K) 

Tc(K) 

820 

293 

331 

900 

293 

341 

820 

318 

358 

900 

318 

369 

820 

343 

377 

900 

343 

391 

Table  5  is  the  test  matrix;  it  presents  a  total  of  nine  (9)  ease  studies  that  have  been  seleeted  for 
validation  purposes  in  this  report.  Cases  A1-A3  eorrespond  to  eonstant  high-pressure  fuel 
injeetion  at  the  standard  Spray  A  eonditions  while  varying  fuel  temperature.  Cases  B1-B3  is  a 
eolleetion  of  intermediate-to-low  vessel  temperature  eondition  at  T^mb  =  820  K  studied  while 
varying  the  injeetion  pressure.  Together  with  eases  B4-B6  they  also  quantify  the  effeet  of 
variations  in  fuel  temperature  at  various  eonditions. 
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Table  5.  CONVERGE  test  cases. 


Case 

Injection  Conditions 

Vessel  Conditions 

CMT-Al 

7>  =  341R;  =  1500  bar, 

Poo  —  60  bar;  Too  —  900K, 
Poo  =  22.8  kg/m'^ 

CMT-A2 

Tf  =  369K-,  Pinj  =  1500  bar 

CMT-A3 

Tf  =  391K-,  Pinj  =  1500  bar 

Case 

Injection  Conditions 

Vessel  Conditions 

CMT-Bl 

Pijij  —  500  bar;  Tf  —  331  K 

Poo  —  60  bar;  Too  —  820^, 
Poo  =  22.8  kg/m'^ 

CMT-B2 

Pi^j  =  1000  bar;  Tf  =  331  K 

CMT-B3 

Pi^j  =  1500  bar;  Tf  =  331  K 

Case 

Injection  Conditions 

Vessel  Conditions 

CMT-B4 

Pinj  —  500  bar;  Tf  =  377  K, 

Poo  —  60  bar;  Too  —  820R, 
Poo  =  22.8  kg/m^ 

CMT-B5 

Pinj  =  1000  bar;  Tf  =  377  K, 

CMT-B6 

Pinj  =  1500  bar;  Tf  =  377  K, 

A  rate-of- injection  boundary  condition  for  the  simulations  was  specified  through  the  use  of 
hydro-dynamic  fuel-mass  flow-rate  application.  This  application  was  developed  by  researchers 
of  CMT-Motores  (Spain)  in  an  effort  to  understand  the  actual  rate  of  injection  during  very  fast 
injector  opening  transients.  It  is  an  online  virtual  flow-rate  generator  that  considers  various 
injector  operating  conditions;  it  is  readily  available  at  http://www.cmt.upv.es/ECN03.aspx.  This 
tool  was  applied  to  generate  the  necessary  rate  of  injections  as  shown  in  figure  7. 


Figure  7.  Rate  of  injection  (ROI)  profiles  for  CMT 

cases  A1-A3  using  the  CMT  virtual  injection 
generator. 
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Figure  8a  shows  that  the  modeled  fuel-temperature  effeets  are  not  as  pronounced  as  the 
experimental  case.  In  measurements,  we  see  about  1-mm  decrease  in  liquid  penetration  length 
when  temperature  is  increased  from  341,  369,  to  391  K.  Simulations  capture  this  trend  but  with 
at  a  reduced  magnitude  of  0.5-mm  per  case.  The  effect  on  the  vapor  penetration  is  negligible  as 
seen  in  figure  8b.  There  were  no  available  measurements  for  comparison  of  figure  8b;  hence, 
modeled  results  are  presented  alone. 


Figure  8.  Spray  penetration  profiles  at  —  60  bar;  —  900K,  —  22.8  kg/vn?  ambient  conditions  with 

fixed  Pijij  —  1500  bar  while  varying  fuel  temperature  7y  =  34:1K,  369K,  391K. 

Figures  9  and  10  show  the  liquid  phase  penetration  response  to  various  operating  conditions  at  a 
fixed  vessel  temperature  ofTamb  —  ^20  K.  It  is  important  to  note  that  this  condition  represents 
the  low-temperature  regime  of  diesel  engine  conditions  and  it  represents  a  departure  from  the 
standard  Spray  A  conditions.  Temperature  variations  in  the  injector  and  vessel  conditions  will 
impact  the  fuel’s  physical  properties  while  imposing  different  heating  loads.  Higher  heating 
loads  and  fuel  temperatures  will  result  in  shorter  steady  state  liquid  penetration  lengths. 
Traditionally  it  has  been  challenging  for  multiphase  CFD  solver  to  capture  low-temperature 
behavior  accurately  and  hence  a  validation  study  in  this  regime  is  important. 

Figure  9a  is  the  rate  of  injection  profile  for  each  injection  pressure  obtained  using  the  CMT 
virtual  injector  generator  for  cases  CMT  B1-B3.  The  total  injected  mass  was  obtained  through 
integration  of  each  profile  using  a  simple  MATLAB  script.  Figure  9b-9d  shows  the  impact  of  the 
injector  pressure  on  the  fuel-penetration  behavior.  As  it  is  expected  higher  rail  pressure  increases 
the  rate  of  atomization  leading  to  modest  decreases  in  liquid  lengths.  Note  that  the  vapor 
penetration  length  increases  with  rail  pressure,  because  it  is  also  dependent  on  the  atomization 
behavior.  In  figure  10  the  same  behavior  is  studied  but  with  a  higher  fuel-temperature  condition, 
cases  CMT  B4-B6.  Figure  10a,  shows  the  rate  of  injection  profiles,  while  figure  lOb-lOd  shows 
the  liquid  and  vapor  penetration  behaviors  with  pressure  at  this  condition.  Because  of  the  steady 
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state  modeling  teehnique  (RANS  turbulenee  model)  no  profile  averaging  was  neeessary,  and  a 
single  shot  study  was  suffieient.  Note  also  that  due  to  the  longer  simulation  time,  4.5  |is  versus 
1.5  gs  for  Sandia  Spray-A  ease,  a  larger  number  of  partieles  were  injeeted  to  avoid  spurious 
oseillations  in  the  profiles.  A  simple  study  determined  that  one  million  total  injeeted  partieles 
were  suffieient  to  aeeount  for  evaporation  and  the  longer  injeetion  proeess.  Also,  the 
eomputational  grid  was  preseribed  based  on  the  resolution  studies  for  Spray  A.  This  was 
determined  using  a  base  grid  eell  size,  dx^ase  —  2  mm,  and  an  embedding  seale  level  of  three,  as 
follows:  dx  =  dxi,ase  ^  ).  Henee  the  minimum  grid  eell  size  at  the  present 

embed  seale  eorresponds  to  dx  =  0.25  mm. 


Figure  9.  Spray  injection  and  penetration  profiles  at  Poo  =  60  bar;  Poo  =  8201^,  Poo  =  22.8  kg/rr?  ambient 

conditions  and  fuel  temperature  Tf  —  331  K.  (a)  ROI  profiles  used  for  cases  b-d,  (b)  liquid/vapor  profiles 
Pin]  —  500  bar,  (c)  liquid/vapor  profiles  at  Pj„y  =  1000  bar,  liquid/vapor  profiles  at  Pj„y  =  1500  bar. 
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Figure  10.  Spray  injection  and  penetration  profiles  at  Poo  =  60  bar;  Poo  =  820K,  Poo  =  22.8  kg/rr?  ambient 

conditions  and  fuel  temperature  Tj  —  377  K.  (a)  ROI  profiles  used  for  cases  b-d,  (b)  liquid/vapor  profiles 
at  Pijij  =  500  bar,  (c)  liquid/vapor  profiles  at  Pj„^-  =  1000  bar,  liquid/vapor  profiles  at  Pj„y  = 

1500  bar. 

In  order  to  enable  comparison  with  current  simulations,  steady  state  results  are  presented  in 
figure  1  la-1  lb  for  cases  CMT  B1-B6.  An  averaged  liquid  penetration  length  was  calculated  in 
the  steady  state  region  from  3.0  to  4.0  ps;  the  same  approach  was  used  with  simulation  results. 
Figure  1  la-1  lb  depicts  the  liquid  length  variations  with  pressure  at  7^  =  331K,  377K, 

respectively.  The  effect  of  fuel  temperature  on  the  penetration  length  is  clearly  observed  at  each 
pressure  interval.  Note  that  the  profiles  retain  their  behavior  while  reducing  almost  uniformly  in 
magnitude.  At  fixed  fuel-temperature  conditions,  the  largest  impact  on  penetration  length  is 
consistently  seen  when  increasing  the  injection  pressure  from  500  to  1000  bar.  Model 
comparison  to  experiments  show  good  agreements  at  each  pressure  interval  for  both  fuel- 
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temperature  conditions.  The  small  discrepancies  between  experiments  and  simulations  may  arise 
due  to  fuel-temperature  calibration  effects.  It  is  also  important  to  note  that  in  experiments  the 
most  accurate  temperature  conditions  was  probed  at  the  injector  sac  volume  as  was  explained 
earlier.  This  temperature  location  reading  is  different  to  what  is  specified  at  the  CFD  boundary 
condition,  which  corresponds  to  the  fuel  temperature  at  the  tip  of  the  nozzle  (downstream  of  the 
sac  volume).  This  discrepancy  in  specification  is  subtle  but  shows  the  need  for  better  consistency 
between  measurements  and  experiments  in  order  to  avoid  calibration  or  tuning  of  CFD  models. 


Figure  11.  Effect  of  7^  variation  over  liquid  penetration  length  for  different  injection  pressures 
(at  =  22.8  kg/m^)  for  (a)  7}  =  331  Tf  ,  (b)  7}  =  377  K. 


5.  Conclusions 


The  studies  conducted  as  a  part  of  this  report  presents  successful  validations  for  nonreacting 
spray  simulations  under  typical  diesel  engine  conditions.  Validation  with  Sandia’s  Spray-A 
unique  operating  condition  was  targeted  to  demonstrate  the  grid-convergence  properties  (P). 

Note  that  AMRs  and  grid-embedding  techniques  were  used  in  obtaining  numerical  convergence. 
Grid-convergence  studies  were  useful  in  determining  an  adequate  grid  resolution  while 
maintaining  model  accuracy.  Subsequent  validation  studies  with  the  CMT-Motores  database 
were  also  targeted  to  explore  the  fuel-temperature  effect  on  spray  parameters  (5).  This  study  was 
conducted  at  a  lower  vessel  temperature  condition  and  represents  a  departure  from  Spray  A 
conditions.  The  effect  of  fuel  temperature  was  seen  to  influence  mainly  the  liquid  penetration 
length  and  to  be  less  pronounced  for  the  vapor  penetration  profiles.  Also,  the  numerical  results 
show  good  agreement  with  measurements.  The  reported  discrepancies  are  currently  thought  to  be 
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due  to  providing  nominal  fuel-temperature  boundary  eondition  in  the  CFD  model  and  the 
“actual”  fuel  temperature  in  the  experiment.  This  finding  raises  the  need  for  better  consistency 
between  measurements  and  experiments  in  order  to  avoid  calibration  or  tuning  of  CFD  models. 

The  results  presented  in  this  briefing  provide  a  benchmarked  model  able  to  predict  engine  spray 
at  various  challenging  diesel  operating  conditions  of  interest  to  the  Army.  Future  efforts  will  be 
aimed  at  conducting  an  extended  range  of  benchmark  studies  with  JP-8  surrogate  fuels  including 
detailed  studies  of  mixture  formation  process,  soot  formation  and  oxidation,  including  split 
injection  strategies. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


AMR 

Adaptive  Mesh  Refinement 

ANL 

Argonne  National  Eaboratory 

ARL 

U.S.  Army  Researeh  Eaboratory 

CFD 

Computational  Eluid  Dynamics 

CPF 

Constant  Pressure  Plow 

CSI 

Convergent  Science,  Inc. 

CT 

computed  tomography 

CVP 

Constant  Volume  Preburn 

DNS 

direct  numerical  simulation 

ECN 

Engine  Combustion  Network 

HEUI 

Hydraulically- actuated  Electronically-controlled  Unit  Injector 

HRJ 

Hydrotreated  Renewable  Jet 

IPK 

Iso-Paraffinic  Koresene 

KH 

Kelvin-Helmholtz 

KH-ACT 

Kelvin-Helmholtz -Aerodynamics  Cavitation  Turbulence 

EES 

large-eddy  simulation 

MPI 

Message  Passing  Interface 

PID 

proportional-integral-derivative 

PISO 

Pressure  Implicit  with  Splitting  of  Operator 

RANS 

Reynolds-averaged  Navier-Stokes 

ROI 

rate  of  injection 

RT 

Rayleigh-Taylor 

UAV 

unmanned  aerial  vehicle 

VTD 

Vehicle  Technology  Directorate 
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